Cosmological dark matter annihilations into 7-rays — a closer look 
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We investigate tlie prospects of detecting weakly interacting massive particle (WIMP) dark matter 
by measuring the contribution to the extragalactic gamma-ray radiation induced, in any dark matter 
halo and at all redshifts, by WIMP pair annihilations into high-energy photons. We perform a 
detailed analysis of the very distinctive spectral features of this signal, recently proposed in a short 
letter by three of the authors: The gamma-ray flux which arises from the decay of vr" mesons 
produced in the fragmentation of annihilation final states shows a severe cutoff close to the value 
of the WIMP mass. An even more spectacular signature appears for the monochromatic gamma- 
ray components, generated by WIMP annihilations into two-body final states containing a photon: 
the combined effect of cosmological redshift and absorption along the line of sight produces sharp 
bumps, peaked at the rest frame energy of the lines and asymmetrically smeared to lower energies. 
The level of the flux depends both on the particle physics scenario for WIMP dark matter (we 
consider, as our template case, the lightest supersymmetric particle in a few supersymmetry breaking 
schemes), and on the question of how dark matter clusters. Uncertainties introduced by the latter 
are thoroughly discussed implementing a realistic model inspired by results of the state-of-the- 
art N-body simulations and semi-analytic modeling in the cold dark matter structure formation 
theory. We also address the question of the potential gamma-ray background originating from active 
galaxies, presenting a novel calculation and critically discussing the assumptions involved and the 
induced uncertainties. Furthermore, we apply a realistic model for the absorption of gamma-rays 
on the optical and near-IR intergalactic radiation field to derive predictions for both the signal and 
background. Comparing the two, we find that there are viable configurations, in the combined 
parameter space defined by the particle physics setup and the structure formation scenario, for 
which the WIMP induced extragalactic gamma-ray signal will be detectable in the new generation 
of gamma-ray telescopes such as GLAST. 

PACS numbers: 95.35.-|-d; 14.80.Ly; 98.70.Rz 



I. INTRODUCTION 



The accumulated evidence for the existence of large amounts of nonbaryonic dark matter in the Universe is by 
now compelling (for a review, see e.g. Data on the cosmic microwave background (CMB)[^ and supernova 

observations ||3[| jointly fix the energy fraction in the form matter and cosmological constant (or something similar) 
to Q.M ~ 0.3 and SIa ^ 0.7, respectively. At the same time, the CMB measurements limit the contribution from 
ordinary baryons to less than U,b ^ 0.05, which is in excellent agreement with big bang nucleosynthesis. This means 
that non-baryonic matter has to make up most of the matter in the Universe, ^dm — ^m- Incidentally, recent 
measurements of the large-scale distribution of galaxies independently confirm JIm = 0.27 ± 0.06 [Q, giving further 
credence to these conclusions. The currently best estimate of Q.m comes from a joint analysis of CMB and large scale 
structure data |^ and gives flMh^ = 0.115 ± 0.009 where h is the Hubble constant in units of 100 km s~^ Mpc~^. 

When it comes to the question of how the dark matter is distributed on small, galactic and sub-galactic, scales 
the situation is much less clear, however (for a review, see, e.g., [^). After being subject to an extensive debate, 
with both theoretical and observational controversies, it seems that the Cold Dark Matter (CDM) model, with dark 
matter made of, e.g., weakly interacting massive particles (WIMPs), or the ACDM model, with a contribution from 
the cosmological constant, are in fair agreement with current observations, so that drastic modifications like strong 
self- interaction are not urgently called for (see, e.g., 0). 

Focusing on the CDM model with WIMPs as dark matter candidates, there is an obvious interest to use as much as 
possible of the knowledge of the distribution of CDM given through state-of-the-art N-body simulations. In particular, 
the distribution of dark matter plays a crucial role in most WIMP detection methods, and determines therefore the 
possibility of identifying the dark matter and pinning down its particle properties. 

In this vein, we recently presented a short note |^ (hereafter BEU) where, contrary to previous predictions it was 
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shown that in the hierarchical picture found in CDM simulations the cosmic 7-ray signal from WIMP annihilations 
may be at the level of current estimates of the extragalactic 7-ray flux. In this paper we deal more carefully with 
the issues of the formation of structure in a CDM or, rather, ACDM universe, investigating the sensitivity of the 
expected gamma-ray flux to different treatments of the structure formation process. We also address the question of 
the diffuse background flux expected from various types of active galaxies and compare its spectral features with those 
of the signal from WIMP annihilations. We consider several sample cases in a theoretically favored WIMP scenario, 
that of supersymmetric dark matter, and highlight the possibility to disentangle such signals from the background in 
future measurements of the the extragalactic 7-ray flux, in particular with the GLAST detector. Results for both the 
signal and background components are presented implementing a careful treatment of the absorption of high energy 
gamma-rays in the intergalactic space caused by pair production on the optical and infrared photon background. 

The paper is organized as follows. In Sec. II we set up the general formalism for computing the dark matter induced 
gamma-ray flux. In Sec. HI we investigate the properties of dark matter halos, on all scales of relevance to our problem, 
in various semi-analytical and numerical simulation scenarios. Implications for the WIMP induced gamma-ray flux 
are discussed in Sec. [V, while in Sec. 0, we investigate the background problem, including the effects of varying 



within present observational limits the slope of the energy spectrum of the gamma-ray emission from active galaxies. 
We also check the effects of removing some more resolved point sources, as may be expected for the next generation 
of experiments. In Sec. Vl we show a few examples of what signals can be expected for one of th e fav ored WIMP 



candidates, the neutralino, and give an estimate of sensitivity curves for the GLAST detector. Sec. VII contains our 
conclusions. 



II. THE DARK MATTER INDUCED EXTRAGALACTIC GAMMA-RAY FLUX 



There are several ways to compute the gamma-ray flux generated in unresolved cosmological dark matter sources. 
In BEU the result was derived by tracing the depletion of dark matter particles with the Boltzmann equation. The 
approach we describe here, in which we simply perform a sum of contributions along a given line of sight (or better, 
a given geodesic), gives the same result but shows more directly the role played by structure in the Universe. We 
assume a standard homogeneous and isotropic cosmology, described by the metric: 



ds^ = c^dt^ - R^{t) [dr^ + Sl{r)dn^] 



(1) 



where dW^ — dO"^ + sin 



and where the function Sk{r) depends on the overall curvature of the Universe: 



r r, fc = 0; 

Sk{r) — < arcsinr, fc = +1; 

I arsinhr, fc = — 1. 



(2) 



In our applications, we will safely use fc = (i.e. we assume a flat geometry for the Universe). The angular interval 
dfl = smdd9d(j) may, e.g., correspond to the angular acceptance of a given detector. At redshift z, dfl and the radial 
increment dr determine the proper volume: 



dV = l^^^^^^drdn 



(3) 



Let dN^/dE[E, M, z) be, on average, the differential energy spectrum for the number of 7-rays emitted, per unit of 
time, in a generic halo of mass M located at redshift z. Even for large M, this source can be safely regarded as 
point-like and unresolved (with the upcoming generation of gamma-ray telescopes, it might be possible to resolve the 
dark-matter induced flux from galaxies in the local group, but almost certainly not further out). Summing over all 
such sources present in dV, we can find the number of photons emitted in this volume and, say, in the time interval 
dt and energy range {E, E + dE); the emission process being isotropic, the corresponding number of photons dN^ 
collected by a detector on earth with effective area dA during the time dto and in the (redshifted) energy range 
(i?o, Ef) + dEo), is equal to: 



dN^ = e-^("'-^°) 



/ dM (M, z) ^ (E, M, z] 

J divr ' ' dE ^ ' ' ' 



dV dA 



47r[i?o5fc(r)]- 



■ dE^ dto 



(4) 



where we applied the relation dtdE = (1 -I- z)^^dto (1 -I- z)dEo — dto dEo, and we introduced the halo mass function 
dn/dM, i.e. the comoving number density of bound objects that have mass M at redshift z (the factor (1 -I- z)"^ 
converts from comoving to proper volume). The first factor in Eq. (|^) is an attenuation factor which accounts for 
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the absorption of 7-rays as they propagate from the source to the detector: the main effect for GeV to TeV 7-rays is 
absorption via pair production on the extragalactic background hght emitted by galaxies in the optical and infrared 
range. Detailed studies of this effect, involving a modeling of galaxy and star formation and a cornparison with data on 
the extragalactic background light, have been performed by several groups (see, e.g., [ij, |l^, |l^ We take 

advantage of the results recently presented by the Santa Cruz group [T5|] ; we implement an analytic parameterization 
of the optical depth t, as a function of both redshift and observed energy, which reproduces within about 10% the 
values for this quantity plotted in Figs. 5 and 7 in Ref. (ACDM model labelled "Kennicut"; the accuracy of 
the parameterization is much better than the spread in the predictions considering alternative models fl^ ). For 
comparison, we have verified that the results presented in Salamon and Stecker jl^ (their model in Fig. 6, with 
metallicity correction) is in fair agreement with the model we are assuming as a reference model in the energy range 
of interest in this work, i.e. below a few hundred GeV. 

The estimate of the diffuse extragalactic 7-ray flux due to the annihilation of dark matter particles is then obtained 
by summing over all contributions in the form in Eq. (Q): 



(IEq dA dVl dto dEo 



^ J dri?oe-^(-^°) J dM^iM,z)^iEoil + z),M,z) 



r I p—r{z,Eo) f j„ 



where the integration along the line of sight has been replaced by one over redshift, Hq is the Hubble parameter, c is 
the speed of light and h depends on the cosmological model. 



h{z) = y/i^Mil + + ^Kil + zf + f^A. (6) 

In this work we put the contribution from curvature — 0, in agreement with the prediction from inflation and 
with recent measurements of the microwave background |^ . Taking the limit in which all structure is erased and dark 
matter is smoothly distributed at all redshifts, Eq. (||) correctly reduces to the analogous formula derived with the 
Boltzmann equation in BEU (Eq. (4) therein). 



III. THE PROPERTIES OF HALOS 



Three ingredients are needed to use Eq. (g) for an actual prediction of the 7-ray flux. We need to specify the WIMP 
pair annihilation cross section and estimate the number of photons emitted per annihilation, as well as the energy 
distribution of these photons: the choice of the particle physics model fixes this element. As photons are emitted in 
the annihilation of two WIMPs, the flux from each source will scale with the square of the WIMP number density in 
the source. The second element needed is then the dark matter density profile in a generic halo of mass M at redshift 
z. Finally we need to know the distribution of sources, i.e. we need an estimate of the halo mass function. 

Some insight on the latter two ingredients comes from the ACDM model for structure formation: we outline here 
hypotheses and results entering the prediction for the dark matter induced flux. We start with the mass function for 
dark matter halos. 



A. The halo mass function 

Press-Schechter theory postulates that the cosmological mass function of dark matter halos can be cast into 
the universal form: 

dn Po dlogv 
dM = Jp''^'^''^dl^ 

where po is the comoving dark matter background density, po ~ Pc^m with pc being the critical density at z = 0. 
We introduced also the parameter v = Ssc{z)/a{M), defined as the ratio between the critical overdensity required for 
collapse in the spherical model Ssc and the quantity a{M), which is the present, linear theory, rms density fluctuation 
in spheres containing a mean mass M. An expression for Sgc is given, e.g., in Ref. jl^. a{AI) is related to the 
fluctuation power spectrum P{k), see e.g. Ref. Q, by: 



a^{M) = / d^k W^{kR)P{k) 



(8) 
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FIG. 1: Fraction of total mass provided by objects heavier than a given mass Al (upper curves) or within 14 decades in mass 
(lower histograms) at three different redshifts and for the mass function as derived in the ellipsoidal collapse model. 

where W is the top-hat window function on the scale = SM/Anp with p the mean (proper) matter density. The 
power spectrum is parametrized as P(fc) oc k"T'^{k); we fix the spectral index n = 1 and take the transfer function T 
as given in the fit by Bardeen et al. iQ for an adiabatic CDM model, with the shape parameter modified to include 
baryonic matter according to the prescription in, e.g. Eq. (15.84) and (15.85). Note that the fit we use agrees 
within 10% with the analytic result obtained for large k in Ref. ||2^, hence holds to the accuracy we are concerned 
about for the small scales we will consider below. We normalize P and a by computing a in spheres of i? = 8//i Mpc 
and setting the result equal to the parameter cg {h is the usual Hubble constant in units of 100 km s~^ Mpc~^). 

In Eq. ([7|) /(j^) is known as the multiplicity function; we implement the form found in the ellipsoidal collapse 
model 

./M^2A(l + -L)(^)"\xp(-if) (9) 

where ly' = y/av, and the parameters q = 0.3 and a — 0.707 are derived by fitting Eq. ^ to the N-body simulation 
results of the Virgo consortium while A is fixed by the requirement that all mass lies in a given halo, i.e. 
/ dh'f{iy) = 1 or J dM Mdn/dM — pQ. Eq. ^ reduces to the form originally proposed in Press-Schechter theory and 
valid for spherical collapse if o = 1, g = and A = 0.5. 

To give the reader a feeling for what the distribution of mass is, as predicted by the halo mass function we are 
considering here, in Fig. |l| we plot the fraction of the total mass in halos heavier than M, and the fraction per mass 
decade, for three different redshifts, z equal to 0, 2 and 4, and for our default choice of cosmological model: = 0.3, 
= 0.7, h = 0.7, Qb = 0.022//i2 and as = 0.73 ||. Note the peak in the distribution at M ~ lO^^ - lO^^ foj. 
z = rapidly moving to lower masses for larger redshifts; note also that the low mass tails are not very steep, with 
only 89% (81%) of the total mass in structures heavier than 10 Af© at z = (z = 2). These numbers get slightly 
larger if one applies the spherical collapse model instead of the ellipsoidal model we have considered here. 

B. The density profile in dark halos 

In the ACDM model for structure formation, dark matter halos are assumed to form hierarchically bottom-up via 
gravitational amplification of initial density fluctuations. Small structures merge into larger and larger halos and final 
configurations are self-similar, with a smooth dark matter component and, possibly, a small fraction of the total mass 
in subhalos which have survived tidal stripping. We neglect for the moment eventual substructure, whose role on the 



5 



WIMP induced signal is discussed in the next Section. N-body simulations seem to indicate that dark matter density 
profiles can be described in the form 

p{r)^p'g{rla), (10) 

where a is a length scale and p' the corresponding density. The function g{x) is found to be more or less universal 
over the whole mass range of the simulated halos, although different functional forms have been claimed in different 
simulations: we will consider the result originally proposed by Navarro, Frenk and White |^ (hereafter NFW profile), 

ffNFW [x] = —-^ — ^2 ' (11) 
x(\ + x) 

supported also by more recent simulations performed by the same group [ p6| , and the result found in the higher 
resolution simulation (but with fewer simulated halos) by Moore et al. (hereafter Moore profile), 

5Moorc {x) = 2,1.5 (1 +2;1.5)- (1^) 

The two functional forms have the same behavior at large radii and they are both singular towards the center of 
the halo, but the Moore profile increases much faster than the NFW profile (non-universal forms, with central cusp 



slopes depending on evolution details have been claimed as well |28 ). There have been a number of reports in the 
literature arg uing that the rotation curves of many small-size disk galaxies rule out divergent dark matter profiles, 
see, e.g., [29[ ^ (note however that this issue is not settled yet, see, e.g., ||3l[), while they can be fitted by profiles 
with a flat density core. We consider then here as a third alternative functional form the Burkert profile p^ , 

g^(^^" (l + x)(l + x^) ' 

which has been shown to be adequate to reproduce a large catalogue of rotation curves of spiral galaxies . 

Rather than by a and p', it is useful to label a dark matter profile by its virial mass M and concentration parameter 
Cvir- For the latter, we adopt here the definition by Bullock et al. |3^: let the virial radius Rmr of a halo of mass M 
at redshift z be defined as the radius within which the mean density of the halo is A„ir times the mean background 
density p{z) at that redshift: 

47r 

M=—A,„p{z)Rl„. (14) 

We take the virial overdensity to be approximated by the expression , valid in a flat cosmology, 

. (187r2 + 82a;-39a;2) 



VIm{z) 

with X = ^Im{z) — 1, (A„ir — 337 for U,m = 0.3 aX z — 0). The concentration parameter is then deflned as 



(15) 



with r_2 the radius at which the effective logarithmic slope of the proflle is —2, i.e. it is the radius set by the equation 
d/dr {r^g{r)) \ ^_^ — 0. This means that r_2 = a for the NFW profile, while X-2 = f_2/a is equal to about 0.63 for 
the Moore profile and to 1.52 for the Burkert profile. Note that these definitions of Rmr and Cmr differ from those 



adopted in Ref. |25| and Ref. [g6|. 

After identifying the behavior in Eq. (^0|), Navarro et al. noticed also that, for a, given cosmology, the halos in 
their simulation at a given redshift show a strong correlation between Cmr and M pq|, with larger concentrations 
in lighter halos. This trend may be intuitively explained by the fact that low-mass halos typically collapsed earlier, 
when the Universe was denser. Bullock et al. [p4l confirmed this behavior with a larger sample of simulated halos 



and propose a toy model to describe it, which improves on the toy model originally outlined in |25 : On average, a 
collapse redshift Zc is assigned to each halo of mass M at the epoch z through the relation Mi,{zc) = FM, where the 
typical collapsing mass is defined implicitly by a {Mi,{z)) = 5sc{z) and is postulated to be a fixed fraction F of M 
(following Ref. |Q we choose F = 0.015). The density of the Universe at Zc is then associated with a characteristic 
density of the halo at z; it follows that, on average, the concentration parameter is given by: 



(17) 
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FIG. 2: Concentration parameter versus mass for halos of mass Al at z — 0. On the left-hand panel we reproduce from Ref. ]34|] 
the behavior found in a large sample of simulated halos, with a binning in mass in which each marker represents the peak in 
the distribution and the relative bar its 68% width; the trend is reproduced with the toy models proposed in Ref. [M itself 
(Bullock et al.) and in Ref. (ENS). On the right-hand side, we show an extrapolation of c„,> to the whole mass range we 
need to include in our analysis according to the two toy models. 



where K is a constant (i.e. independent of M and cosmology) to be fitted to the results of the simulations. Bullock 
et al. show that this toy model reproduces rather accurately the dependence of c^r found in the simulations 
on both M and z. We reproduce this fit at z = in Fig. ^ (left panel, solid line); "data" points and relative error 
bars are taken from |34[ and just represent a binning in mass of results in their simulated halos: in each mass bin, 
the marker and the error bars correspond, respectively, to the peak and the 68% width in the Cyir distribution. We 
determine K with a best fitting procedure in the cosmology fin = 0.3, f^A = 0.7, h = 0.7 and as = 1 adopted in 
the N-body simulation referred to, and then use this value to estimate the mean Cmr in other cosmologies; we find 
K = 4.4. Finally, following again Bullock et al. |Q, we assume that, for a given M, the distribution of concentration 
parameters V is log-normal with a Icr deviation A(log]^o Cuir) around the mean, independent of M and cosmology; we 



take A(logj 



.) =0.2. 



An alternative toy-model to describe the relation between Cyir and M has been discussed by Eke, Navarro and 
Steinmetz (hereafter ENS model): The relation they propose has a similar scaling in z, with however a different 
definition of the collapse redshift Zc and a milder dependence of Cyir on M. In our notation, they define of Zc through 
the equation: 



L»(zc)CTcff(Mp) 



1 

c: 



(18) 



where D{z) represents the linear theory growth factor, and acs is an 'effective' amplitude of the power spectrum on 
scale M: 



dln(CT) 



(M) 



da 



-M 



(19) 



dln(M) 

which modulates a{M) and makes dependent on both the amplitude and the shape of the power spectrum, rather 
than just on the amplitude as in the toy model of Bullock et al. Finally, in Eq. (p^, Mp is assumed to be the mass of 
the halo contained within the radius at which the circular velocity reaches its maximum, while Ca is the parameter 
(independent on M and cosmology) which has to be fitted to the simulations. With this definition of Zc it follows 
that, on average, Cvir can be be expressed as: 

1/3 

Cy„{M,z)^\- — . (20) 



.(z) U.m{Zc) 



1 + z 
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As we already mentioned the dependence of Cmr on M as given in the equation above is weaker than in the Bullock 
et al. toy- model. Our best fitting procedure gives = 76 and the behavior in Fig. || (left panel, dashed line), which 
reproduces the N-body "data" fairly well, with values not very far from those obtained in the Bullock et al. model 
within the range of simulated masses, and possibly just a slight underestimate of the mean value in the lighter mass 
end. 

On the other hand, the extrapolation outside the simulated mass range can give much larger discrepancies as shown 
in the right panel of Fig. ^. Solid lines are for the same models as those shown in the left panel {K and Co- from 
the data fit in the left panel), with just erg set equal to our preferred value, as — 0.73. When going to small M , Cmr 
increases in both cases, but the growth in the model of Bullock et al. is much faster than in the ENS model; we 
will show explicitly how this uncertainty propagates to the prediction of the dark matter induced 7-ray flux. The 
sensitivity of our results to the choice of cosmological parameters is generally much weaker: The largest effect is given 
by the overall linear scaling of Cyir{M, z) with erg- There is also the possibility to change the cosmological model by 
including other dark components; we are not going to discuss any such case in detail, we just mention that a neutrino 
component at the level of current upper limits is not going to change severely our picture, while a substantial warm 
dark matter component may play a crucial role if Zc is indeed defined according to the ENS prescription. 



IV. WIMP INDUCED FLUX: ROLE OF STRUCTURES AND SPECTRAL FEATURES 



We are now ready to write explicitly the term dAf^/dE introduced above and to derive the formula for the flux. 
The differential energy spectrum for the number of photons emitted inside a halo with mass M at redshift z is 



'dE 



{E,M,z) 



av dN^jE) 
"2 dE 
avdN^{E) M 



dc' 



A,; 



dE 



Ml 



dc' 



d^r 



\r/a) 



X-2) 



■I2{X. 



ram 7 "-^vir 



(21) 



where av is the WIMP annihilation rate, dN^{E) / dE is the differential gamma-ray yield per annihilation and 
is the WIMP mass. Note that in previous literature, the prefactor av/2 has often been erroneously taken as av. 
The derivation based on the Boltzmann equation in BEU adds to our confidence that the factor should he av/2. In 
Eq. (pi|) we applied the definition of Rvir and introduced the integrals 



dxx^g'\x) 



(22) 



with the lower limit of integration Xmin = fmin/o'' set, in a singular halo profile, by WIMP self annihilations, i.e. 
roughly by p{rmin) — rn^/[av (to — tc)], where to is the age of the Universe and tc is the collapse time for the halo 
under investigation. To include all sources labeled by their mass M, we averaged over the log-normal distribution 
V{clir) centered on Cvir as given in Eq. ( |l7| ) or (|20|). 

Inserting Eq. ( ^l| ) into Eq. (|^) , we find that the gamma-ray flux is 



dEo 



av c Pq 



dz (i^,)3Ajz)diV.(i^o(l + .))^_.(,.„) 
h[zj dE 



where we have defined 



and the quantity 



A2 



M 



{z,M) 



A'iz) 



dM 



y{z,M)f{v{z,M)) 
a{M) 



da 



dM 



dc' Vic' ^ 



hiXmin, C^^ir X-2) 



(Z,M)X_2)]' 



(c^,(z,M)x_2)' 



(23) 



(24) 



(25) 



(Note that this definition differs from that in BEU js) by a factor 1/(1 -I- z)^. The advantage of the present definition 
is that A\.j{z,M) = 1 if all matter is at the mean density for redshift z.) In early estimates of the WIMP induced 
extragalactic 7-ray flux, see, e.g., the role of structure was not appreciated and the dark matter distribution was 
assumed to be described simply by the mean cosmological matter density p{z) = Pc^m{^ + z)^- Compared to this 
picture, A^j(z, M) gives the average enhancement in the flux due to a halo of mass M, while A^(z) is the sum over 
all such contributions weighted over the mass function. As we will see, the enhancement of the annihilation rate due 
to structure amounts to several orders of magnitude. 
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FIG. 3; Average enhancement in the 7-ray flux emitted in a halo of mass M at redshift z — with respect to the case in which 
the same amount of dark matter is smoothly distributed. On the left-hand side we show how sensitive the result is to the 
concentration parameter. On the right-hand side the result for three different families of dark matter density profiles is shown. 



A. Flux normalization 



We analyze first how sensitive the flux is to the dark halo properties we discussed in Section [11 . In Fig. ^ - left 
panel - we plot A|j as a function of M, at redshift z = and assuming the Moore density profile to describe dark 
matter halos. The four cases displayed correspond to the two toy models for Cmr we have discussed in the previous 
Section, and to two choices of as- our default value as = 0.73 |Q and the larger value as = 1.02 found in another 
recent analysis |3|] and more on the line with values often quoted in the past. For small M, i.e. large Cmr, A\j scales 
roughly like c^j^/ log (cmr), where the logarithmic term follows from the fact that the halo profiles we are considering 
have logarithmically divergent masses (which we cut at the virial radius). It follows that the uncertainty on erg induces 
about a factor of 2 uncertainty on , while an indetermination of a factor of a few is due to the model applied to 
extrapolate Cvir to small masses. 

In Fig. ^ - right panel - we restrict to Cmr as computed with the Bullock et al. toy model, and show the dependence 
of the signal on the choice of halo profile. The spread in the predictions between the Moore profile and the Burkert 
profile is around a factor of ten independent of mass, which is much smaller than the uncertainty due to the choice of 
profile when considering the dark matter induced 7-ray flux generated in single resolved sources. This is one of the 
advantages of considering the cosmological signal. Of course, one is also less sensitive to the actual halo properties of 
a single galaxy, the Milky Way, which are poorly known. This issue is analyzed further in Fig. ^ where, for a given 
halo of density profile p(f), we plot the dimensionlcss quantity 



AS2 8.5kpc VO-SGeV/cm^y J^^ J^^^ 

a sum over contributions along the line of sight in a cone of aperture Afl in the direction ^ (this quantity often appears 
in analyses of the WIMP induced flux generated in the Milky Way halo; normalization factors are fixed following the 
choice in Ref. |^). We focus on a 10^^ Mq halo, i.e. a halo of the size of the Milky Way or Andromeda, and assign 
to it the mean Cyir in the Bullock et al model; also, we choose ip in the direction of the center of the halo and consider 
a moderately large acceptance angle. Ail = 10^'^ sr. We let then the distance d between the center of the halo and 
us vary between 10~^ R^r and 10^ Rmr [Rvir ^ 260 kpc in our sample case) and plot the corresponding (J) for the 
three halo profiles we introduced. The arrows in the figure marks the location on the horizontal axis of the Milky Way 
(MW) and Andromeda (M31). At large d/Rmr we find for all halo profiles the scaling one expects for point-like 
sources: this is obvious for ratios larger than d/R^r ~ \/ tt/AH, ~ 56, when the halo is fully contained in the field of 
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FIG. 4; Scaling of the collected 7-ray flux with the distance d between the detector and the center of a halo, for three difl^erent 
halo profiles. The angular acceptance of the detector is assumed to be Af2 — lO'^^sr. The plot is for a lO^^M© halo, the 
arrows indicate the position on the horizontal axis for the Milky Way and Andromeda; the case for other masses is analogous. 



view; however, as it can be seen, for the Burkert and the NFW profiles such scaling appears already for ratios one 
order of magnitude smaller, and it is present essentially over the whole range displayed for the Moore profile. This 
indicates that the bulk of the flux is emitted in the inner halos: for the Moore profile 50% (10%) of the total emitted 
flux is generated within a radius that is about 9 • lQ~^Rvir (6 • 10~^i?„ir), for the NFW and Burkert profiles the 
corresponding radii are shifted, respectively, to 2.4 • 10~^i?„ir (3.3 • 10~^i?^ir) and 6 • lQ~'^Ryir (2.4 • 10~^i?„ir). While 
the spread in predictions for the flux generated in the center of our Galaxy is very large (6 orders of magnitude), the 
total emitted flux is a much weaker function of the density profile - the uncertainty is roughly an order of magnitude. 

This factor of 10 uncertainty is nearly independent of M , therefore it propagates as an order of magnitude uncer- 
tainty on the overall normalization of the WIMP induced 7-ray flux. 

The behavior of is obtained by folding the scaling of the integrated mass function in Fig. |l| with that of 
in Fig. 1^. The dominant contribution to A^ comes from very small halos: the integrand in A^ is the product of two 
mildly divergent quantities, the mass function times M and A^^; the result is still convergent but heavily relies on 
our understanding of the light mass end. This is shown in the left-hand panel of Fig. ||, where, for the Moore profile 
and our preferred cosmology, we plot A^ at z = restricting the range of integration over mass. For the Bullock et al. 
toy model the contribution per logarithmic interval keeps increasing even for the lightest mass range displayed, while 
in the ENS model it starts decreasing but rather slowly. Extrapolations of Cmr with our toy models to exceedingly 
small masses may not be fully reliable; we prefer to introduce a cutoff in Cyir and hence in A^ at the intermediate 
mass range Mcut, say IO^Mq for z = 0, where we believe the toy models are sufficiently trustworthy. We assume: 

Cy^r{M,z)=C„r{Mcut,z) \f M < M^ut (27) 

The choice of Mcut is to some extent arbitrary; should one make a different assumption Fig. |5| allows to scale our final 
results. 

In Fig. - right-hand panel, we plot (1 -I- z)^A'^/h, i.e. the quantity we need to integrate over z to get the 7-ray flux, 
see Eq. (p3), once folded with the emission spectrum and the absorption factor. We consider both models for Cmr 
and two schemes to define Mcut- In the first we fix Aleut — W^Mq for any z, progressively reducing the mass range 
over which Cyir is extrapolated. Another possibility is to keep the range of this extrapolation fixed: at z—0 we choose 
Mi,/Mcut^ with Aii, the largest scale allowed defined implicitly by a {Mi,{z)) = 5sc{z) and again Aleut = 10^ M© ; at 
other z the same ratio is imposed (we never include extrapolations to masses lower than lOM© ; at the redshift of a 
few when Mcut would be lower than that, we set Mcut = lOAf© ). Both schemes are rather arbitrary, we will show 
however that the final result is not very sensitive to them. Notice, on the other hand, the sharp increase of (1 -I- 2)'^A^ 
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FIG. 5: Enhancement in the diffuse 7-ray flux compared to the case when all structures in the Universe are erased. On the 
left-hand side the contributions of structures of given masses at z = are shown; on the right-hand panel we show the redshift 
dependence, rescaled with the term (1 -I- z)"^ /h{z). 



at small z for the ENS model, whereas a mild decrease or a flat behavior is fomid in the Bullock et al. model. At 
larger z, the scaling in l/h{z) takes over. 



B. Spectral signatures 

We now try to estimate, in an approximate way, the level and spectral shape of the gamma-ray flux that can be 
expected for a general WIMP, leaving a more detailed discussion of the extragalactic background one has to fight 
against for Section 0, and predicted signals for a more specific (supcrsymmetric) dark matter candidate for Section Vl . 

The differential gamma-ray yield per WIMP pair annihilation can be written as: 

= ^ b,^n,x5 {E - (1 - Mi/4M^)) + ^ h/-^{E) . (28) 

X F 

The first term refers to prompt annihilation into two-body final states containing a photon, which, forbidden at 
tree-level essentially by definition of dark matter (zero electric charge), are allowed at higher order in perturbation 
theory. Although subdominant, they have the peculiarity of giving monochromatic 7-rays: as WIMPs in halos are 
non-relativistic the energy of the outgoing photon is fixed by the WIMP mass and the mass of the particle X (i.e., 
E = for the 27 final state and E — M^{1 — M^/4M^) for final states with some non-zero mass particle X). The 
parameter b-yx is the branching ratio into these channels and n-yx is the number of photons per annihilation, i.e. 2 for 



the 27 final state and to 1 for the others. The second term in Eq. (28) is instead the term due to WIMP annihilations 
into the full set of tree-level final states F, containing fermions, gauge or Higgs bosons, whose fragmentation/decay 
chain generates photons; this process gives rise to a continuous energy spectrum. 

Although there is some span in the predictions for the photon emission rate in different particle physics models, 
the spectral features of the induced fluxes are quite generic and can be outlined without referring to a specific model 



(in section VI below we will discuss results for more specific models). We start discussing the monochromatic terms, 
focusing to be definite on the process XX ~^ ^7 and picking for reference some typical value for the annihilation cross 
section in this channel. Consider, e.g., that in the simplest case (no resonances or thresholds near the kinematically 
released energy in the annihilation 2M^) the WIMP total annihilation rate is fixed by the approximate relation [[40[ : 

av ^ (av) ^ 3 ■ 10-^^ Ws'^ ^ ^ iQ-^^^cmh-' , (29) 
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FIG. 6: Spectral signature in the extragalactic gamma-ray flux due to the annihilation of dark matter WIMPs into monochro- 
matic photons. A toy model with three choices of WIMP masses, M-^ = 50, 100, 250 GeV, and fixed annihilation rate into 27, 
is considered. The signature arises because of the asymmetric distortion of the line due to the cosmological redshift, as well as, 
by absorption of gamma-rays generated in distant sources. The normalization of the fluxes are computed assuming halos are 
modelled by the Moore profile and concentration parameters are derived with the Bullock et al. toy model (left panel) or the 
ENS model (right panel); solid and dashed curves refers to two schemes for the choice of the halo mass cutoff Mcut, which, as 
can be seen, plays a marginal role. 



which shows the order of magnitude scaling between the thermaUy averaged annihilation cross section (av) and the 
WIMP thermal relic abundance O^. Note that this relation is only a rough approximation and that large deviations 
from it can appear mainly due to resonances and thresholds. In Section VI below we will not use this approximate 
relation, but instead calculate the relic density including properly both coannihilations, resonances and thresholds. 
For the current discussion though, this approximate relation suffices. The annihilation into two photons is a 1-loop 
process so, in general, its strength is much smaller than av; we assume, as a sample case when this channel is relevant, 
62^ = 10-3. 

In Fig. ^ we show the induced extragalactic gamma ray flux for 3 different values of the WIMP mass, = 
50, 100,250 GeV, and for the two schemes we have considered to estimate Cmr- We c onside r halos modeled by the 
Moore profile, with no subhalos (the effects of the latter will be discussed in Section IV C). The figure illustrates 
the novel signature, first proposed in BEU, to identify a WIMP induced component in the measured extragalactic 
gamma-ray background, the sudden drop of the gamma-ray intensity at an energy corresponding to the WIMP mass 
due to the asymmetric distortion of the line caused by the cosmological redshift. The energy of the 7-rays at emission 
determines whether the smearing to lower energies has a sharper or smoother cutoff: for a larger the absorption 
on the extragalactic optical and infrared starlight background becomes more efficient. Spectra obtained applying the 
ENS toy-model for Cmr are similar to those derived with the Bullock et al. model; the main difference, for masses 
lower than about 100 GeV, is a slight shift of the fiux peak to lower energies. This effect, due to the sharp increase in 
shown in Fig. ^ in a regime where the absorption factor does not rapidly take over, tends to reduce the difference 
in the fiux normalization one might have foreseen looking at alone (in the next generation of measurements the 
energy resolution will probably not be better than 10% or so). Fig. || also illustrates the fact that, at least for the 
line contributions, the treatment of the cutoff in halo mass is not very important; there is mainly an overall scaling 
with the choice of Mcut at z = 0, which the reader can infer from the left-hand panel of Fig. ^ 

Typical features in the continuum contribution are illustrated in Fig. |7[ We have assumed again that av ^ 
3 • lO-^^cm^ and supposed that, as often happens in real particle physics models, the dominant annihilation 
channel is into bb; the energy distribution per annihilating pair in their rest frame is simulated with the Pythia 
Monte Carlo package As most photons are produced in the hadronization and decay of 7r"s (98.8% decay mode: 
7r° — > 27), the shape of the photon spectrum is peaked at half the mass of the pion, about 70 MeV, and is symmetric 
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FIG. 7: Spectral features for the extragalactic gamma-ray flux due to the photons with continuum energy spectrum emitted 
in pair annihilations of dark matter WIMPs. The cases considered here are for a WIMP toy model of given mass and flxed 
total annihilation rate, assuming the dominant branching ratio is bb. The flux normalizations are computed under the same 
assumptions as in Fig. ^. In the left panel we compare the shape of the induced flux at Earth with the one at emission; in 
the right panel we show the dependence of the spectral shape on WIMP mass. For comparison, the EGRET estimate of the 
extragalactic background flux is shown. 



around it on a logarithmic scale (sometime this feature is called the "tt^ bump", see, e.g., |4^). Had we chosen 
a different dominant annihilation channel or a combination of several channels, we would have found very similar 
behaviors. As absorption becomes negligible going to low energies, two features arise when summing extragalactic 
contributions over all redshifts: The peak in the spectrum is shifted to lower energies and there is a sharper decrease 
in the flux approaching the value of the WIMP mass. The first signature is probably hidden in background fluxes, 
see the discussion in the next Section. The second feature is instead potentially interesting, especially in case the line 
components are negligible: While a sensible contribution to the extragalactic 7-ray background can be provided by 
WIMPs in the few GeV energy range, at higher energies the WIMP induced fiux is very rapidly suppressed. Such 
behavior cannot be associated to a spectral index, while background components are closer to a power law. 

As shown on the right-hand side of Fig. |^, the WIMP-induced extragalactic flux gradually flattens for heavier and 
heavier WIMPs; also shown is the current estimate of the diffuse extragalactic background flux as derived from the 
analysis of data taken by the EGRET telescope p3|] . 



C. Role of subhalos 



We have shown that small dense halos are providing the bulk of the WIMP induced 7-ray flux. So far we have 
considered just the case for isolated halos; as already mentioned, N-body simulations indicate that, in the clustering 
process with large halos forming by the merging of smaller objects, a fraction of the latter, up to about 10% of the 
total mass, may have survived tidal disruption and appear as bound subhalos inside virialized halos From 
the point of view of structure formation, the presence of rich substructure populations was at first seen as the main 
flaw in the picture from N-body simulations of ACDM cosmologies, a "crisis" urging for a solution maybe 
with a drastic change in the particle physics set up, see, e.g., Q. More recent analyses indicate that those results 
should be reinterpreted and the apparent discrepancy between the number of subhalos found in the simulations and 



that of luminous satellites observed in real galaxies is fading away, see, e.g., |4^, 49, From the point of view 
of dark matter detection, substructure may play a crucial role ||51| , even in the interpretation of currently available 
data. Consider, e.g., the gamma-ray halo surrounding the Galaxy for which statistical evidence has been claimed in 
data collected by the EGRET telescope the conjecture that this may be generated by pair annihilations of relic 
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dark matter particles is based on the possible presence of dense dark matter clumps in the Milky Way halo | |53| . 

At first sight, the role of substructure may seem marginal in our context. The fraction of mass in subhalos is 
small and the subhalo mass function is not likely to be significantly steeper than the mass function for isolated halos, 
Eq. @). We expect then the number of halos in a given mass range to be larger than the number of substructures 
in the same range. On the other hand, the concentration parameter in subhalos may be significantly larger than for 
halos: on average, subhalos arise in higher density environments, as well as we expect a depletion in their outskirts 
by tidal stripping. This trend has indeed been observed in the numerical simulation of Ref. [Q, where it is shown 
that, on average and for M ~ 5 • IO^^Mq objects, the concentration parameter in subhalos is about a factor of 1.5 
larger than for halos. 

Consider a halo of mass M and suppose that, on average, a fraction / of its total mass is provided by substructures 
with mass function dng/dMg. The differential energy spectrum for the number of photons emitted in such a halo, 
rather than by Eq. (pl|), is now given by: 



dE 



dE 



dc' Vic' ) 



d^r g {r/a) 



.M.i^/.c;.,.P.,o;.,.)(|-)7.'--/ 



av dNj (E) pc 
~2 dE m| 



rg {r/a) 



dn^ 



(1 - f)^MAljiz, M) + I dM.^MsAli^ {z, Ms) 



(30) 



A simple ansatz is that the subhalo mass function has a power-law behavior dus/dAIg oc 1/Aff for Mg < M {(3 < 2 
is required for the total mass to be finite), with the normalization fixed by using the definition of /, i.e. 



(31) 



This gives 



(32) 



If we further assume that / and dug/dMs are independent of M, inserting Eq. ( p2| ) 
contribution of subhalos can be included in the formula for the gamma-ray flux, Eq. 



into Eq. (|3^), we find that the 
31), with the replacement: 



Al,(z,M) ^ (l-/)2 Ai,(^,M)H-(2-/3)/M'3-2 / dM^Ml'^ Al,^{z,Ms) 



(33) 



Here /S.\j is just but with values of Cvir and 7'(c^j^) appropriate for the subhalos. It may be premature to 
deduce the latter from N-body simulation results. The scaling Cmr oc M~'^-'^ proposed in Ref. Q probably cannot 
be extrapolated to small masses: for lO^M© subhalos we would get a value of the concentration parameter 40 times 
larger than the value for halos as computed with the Bullock et al. toy model. 

A prediction for the subhalo mass function is missing as well; there are just limited studies, not fully covering the 
mass range we are interested in. The current N-body simulation results are consistent with a power law behavior, but 
with non-universal slope and some indication that the index (3 is getting harder decreasing the mass of the host halo. 
We find, e.g., from Fig. 5 in Ref. |Q that B ^ 1.66 for a lO^^M© halo. A few studies are focused on Milky Way size 
halos, lO^^Af© : from, e.g.. Fig. 1 in Ref. |Q we can extract the scaling dng/dM^ oc l/M^-^'^. 

The value of / is a matter of debate as well. Values for the fraction of mass in substructures quoted in the literature 
are in the range 1 % to 10 % and they often refer to the ratio of the sum of the masses of identified subhalos to the 
total mass, rather than to an extrapolation performed assuming a mass function. Such a value, say /', should then 
depend both on the algorithm for finding subhalos in the simulation, and, most importantly, on the resolution of the 
simulation. Suppose that /' refers to a simulation where, for halos of mass M, substructures of mass down to pAI 
can be resolved; then, with our notation. 



f'M = /Af^-i [M 



2-13 



(pM) 



2-/31 



(34) 



i.e., / = /7(1 If M = IO^^Mq and p = 5 • 10"^ g|, we find / = 2.56/'. It is then not implausible that the 

true /, eventually to be found at future ultra- high resolution simulations, may approach or even exceed 10 %. 
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FIG. 8; Influence of substructure on the flux normalization for three different average fractions / of the total mass in subhalos; 
we have restricted to a speciflc mass function (see text) with spectral index 13 and kept as a free parameter, we display on the 
horizontal axis the mean enhancement in the concentration parameter in subhalos. 



To give a feeling for the possible effect of substructure, we consider the simplified sample case in which / and the 
mass function are universal, and keep the average enhancement in the concentration parameter as a free parameter 
(we find the value of c„ir(M, z) for subhalos by a rigid rescaling of the Cyir{M, z) found for halos of same mass and at 
the same z: actually, the mass range in which this rescaling matters is just around the cutoff mass Mcut = IO^Mq ). 
In Fig. H we consider (3 = 1.95 or the slightly softer (3 = 1.90, choose three sample values for the fraction of the 
mass in subhalos / and plot the ratio of the value of with and without including subhalos as a function of the 
average enhancement in the concentration parameter. Sensible gains in and hence in the 7-ray flux normalization 
are viable even for moderate enhancements in the concentration parameter. Again, the effect of substructure is less 
dramatic than in case of single dark matter sources: the argument here is analogous to the one presented in the 
discussion on the role of the singularity in halo profiles. 



D. Observability of subhalos in the Milky Way halo 

It would be of utmost importance to test the subhalo picture predicted by CDM N-body simulations by collecting 
information from the morphology of the Milky Way halo. As already mentioned, a rich population of luminous 
satellites is not observed in the Galaxy and this was considered, up to recent work, one the most severe "problems" 
of CDM. There are now models [^9[ |5^ to explain why small substructures may be totally dark (without visible 
baryons); if this is indeed the case, WIMP annihilation might be the only chance to perform a detailed mapping 
of the distribution of mass in the Milky Way. This issue has been investigated by numerous authors (for a recent 
analysis see, e.g., Ref. Q). The problem however reduces to the study of the actual realization of incalculable random 
processes and this implies that it is very hard to estimate the probability for detection of a signal. In particular, a 
crucial parameter will be the location of the nearest dark matter clump, since this will dominate the signal. We show 
here how the picture we have outlined for a generic halo applies to the Milky Way and discuss the implications for 
the observability of subhalos. 

The gamma-ray flux from a single "clump" of mass Ms and at the distance d from the Earth is equal to: 

dE 2 dE ATTj^n Ji.o.s. \M^J 
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FIG. 9: Maximum distance dMax of a clump from the our position in the Galaxy for which the 7-ray flux due to WIMP 
annihilations in the clump exceeds the estimated point source sensitivity of GLAST. We picked a specific WIMP dark matter 
candidate, while we are considering a few models to relate the mass of the clump Ms to its concentration parameter, as well as 
two models for the halo density profiles in the clump (see the text for details) 



The angular extension of most clumps is much smaUer than Afi, hence we can use the point source approximation. 
In our formalism the formula becomes: 

_ (Jv dN^jE) 1 Ms KrrP^ {cUrX-2f , . ,^r,. 

dE - 2 dE ^.d^Ml 3 ^''^ 

Assume (\>^f{E > 100 McV) > 1.6 • lO^^cm^^s"^ is the point source sensitivity of the next gamma-ray telescope 
in space, GLAST, which, as EGRET did, will map the whole gamma-ray sky. In defining a particle physics model 
for a WIMP, one has to fix M^, av and the branching ratios into each annihilation channel. It is then possible to 
compute, say, 

J 100 McV "-f^ 

For each Mg we can estimate the maximum distance d]^Iax of the clump from us such that the WIMP induced flux is 
larger than the point source sensitivity of GLAST. This is shown in Fig. ^ for one of the WIMP toy models introduced 



in Ref. fVB: we assume that = 50GeV and that the total annihilation rate into bb is av — 10 cm'^s , and 



find iV^'-"' = 25.9; a generalization to other models can be obtained very simply by scaling of these values. Six 
configurations for the normalization of the flux are chosen: we assume that Cvir {M, z = 0) for subhalos is either equal 
to the mean value found with the Bullock et al. or ENS toy models for isolated halos (labels "B. et al." and "ENS" 
respectively) or to 4 times the value found with the Bullock et al. model (label "subh."); we consider also the cases 
that subhalos are described both by the Moore et al. profile and by the NEW profile (the results for the Burkert profile 
are very close to the ones for NEW, unless the clump is very close to us, see Eig. ||). As can be seen, unfortunately, the 
prediction of our analysis is that just a few nearby clumps might be detected by GLAST. Eor comparison, we show 
the location in the plane of the figure of a "clump" that is sufficiently massive to have a luminous counterpart, Draco. 
This dwarf spheroidal, together with other similar candidates, has been considered several times in the literature as 
a potential gamma-ray dark matter source, see, e.g., [57| (note, however, that our picture applies on average, rather 
than to a single specific source, which might be better characterized through its rotation curve). 
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V. THE DIFFUSE EXTRAGALACTIC GAMMA-RAY BACKGROUND 



The observability of the signal proposed here depends on the level of the diffuse extragalactic gamma ray background. 
Contributions from several classes of unresolved discrete sources have been discussed in the literature. After EGRET 
maps of the gamma-ray sky, the case for a dominant contribution from blazars is generally considered to be very 
strong: a large number of high-energy emitting blazars has been observed and, as we will show and contrary to other 
candidates, their distribution of energy spectra seems to be compatible with the observed extragalactic radiation. We 
will then rederive here the expected diffuse background under the assumption that the source of the background is 
unresolved blazars. We will mostly follow the analysis of Salamon and Stecker |^, |5^, |6^, but update it with more 
recent data and examine the expected uncertainties. 



A. Basic blazar model 



The basic model assumes that the diffuse gamma ray flux comes from unresolved blazars. We will assume that 
the blazars are distributed in redshift and luminosity according to a luminosity function (Pj , z) where Pj is the 
luminosity (in units of W ]iz~^ sr~^). The luminosity function is the comoving density in units of Mpc~'^ (unit 
interval of logj^Q P^/^)"^. We will further assume that the blazars emit gamma rays with some spectral index a which 
is distributed according to a distribution function p{a). The absorption of gamma rays emitted at redshift z and 
observed at energy Eq is, as before, parameterized in terms of the optical depth t(Eo, z) such that the attenuation 
is proportional to e~'^^^°'^\ We will here use the parameterization of the Kennicut model in Primack et al. |l^ 
introduced in Section II. For comparison we will also use the model of Salamon and Stecker ( their Fig. 6 with 
metallicity correction) . There is also a recent estimate of the absorption by de Jager and Stecker p^ , but we will not 
use that model since it is not valid above z ~ 0.3 which is not sufficient for our purpose. In EGRET observations one 
has seen that blazars most of the time are in a quiescent state but some small fraction of time are in a flaring state 
with higher luminosity and slightly different spectral index (softer, i.e. higher a). We will assume that the blazars are 
in the flaring state a fraction C of their time and that their luminosity then is a factor A f higher than in the quiescent 
state. We will also parameterize the different spectral indices by assuming that they come from the same distribution 
function p{a) but shifted by Aa^ and Aa/ for the quiescent and flaring states respectively (these two quantities are 
not independent; we have adopted here the same notation as in Salamon and Stecker, but, alternatively, one could 
redefine a and introduce a single shift Aa). Putting this together we can write the gamma ray flux (in units of cm^^ 
sr~^ GeV""'^) at energy Eq as 



^7(^0) = jr dlOg^gP-y 



da 



dz 
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P^ 



h{z) 2'KnE^j 
In this equation, we have introduced the following 
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dE 
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= Hubble constant today 

= cosmology factor as defined in Eq. (6) 

0.1 GeV ~ the fiducial gamma ray energy at which the luminosity is P-y 

= the gamma ray spectrum (normalized to 1 at E^j) 

= Planck's reduced constant 

= the luminosity function for blazars in quiescent state 

= the luminosity function for blazars in flaring state 



(39) 



Note that we have for clarity explicitly included c and fi in Eq. (p8|), but the unit conversion factors to get the flux 
in the above given units are not given explicitly. Note that there is no factor of l/47r since the luminosity P^ is per 
sr already. We will, as before, assume that Hq = 70 km s~^ Mpc~^, Vim = 0.3 and Oa = 0.7. 
In the following subsections we will go through the different factors entering in Eq. (|3J 



B. The luminosity function p^{P~,,z) 

We need to know the luminosity function as a function of redshift. Since not that many blazars are observed we 
will follow pq, p9|, pO| and assume that the same basic mechanism (i.e. the same population of high-energy electrons) 
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is responsible for both the gamma ray and the radio flux. We can then use the much larger catalogs of radio sources 
to get the luminosity function. We will assume that the luminosities in gamma and radio are related by 



P. 



P. 



7,/ 



AfKPr 



(40) 



where P.y and P,. are the luminosities (in units of W Hz^^ sr^^). The gamma ray luminosity is given at 0.1 GeV and 
the radio luminosity at 2.7 GHz. The subscripts q and / refer to the quiescent and flaring states respectively. We will 
assume that the two luminosity functions are related by 



P-f{P-f, z) = r]pr{Pr, z) 



(41) 



where p.y and pr are the luminosity functions (in units of Mpc^'^ (unit interval of logj^g P-y/r)^^)- The factor rj takes 
into account possible beaming effects which could mean that not all radio blazars emit gamma rays towards the Earth 
(or vice versa). Including the effect that the blazars are assumed to be in the flaring state a fraction C, of the time, 
we can finally write 
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For the radio luminosity function, we use the parameterization by Dunlop and Peacock [ ]6l[ 
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(43) 



valid up to 2; = 5. This luminosity function was derived for a cosmology with Hq = 50 km s^^ Mpc~^ and ilo = 
f2m = 1, but we can approximately convert this to a luminosity function for our cosmology by multiplying p with a 
correction factor Ipll 



Hoh(z) 



dV 



and multiplying Pdz) with the correction factor 
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^ ^ / our cosmology 



(44) 



our cosmology 



{{1 + z)RoSk{r)) 

our cosmology 
((1 + ^)^0'5'/c(r))^^^^^^;L,//o-50 kms-i Mpc-i 



(45) 



where is the luminosity distance. The luminosity function Eq. (Esf) is valid between P"™'" = 10^® W Hz ^ sr ^ 



and P^^^^^ — 10^° W Hz~^ sr~^ which we will convert to limits on Pj. Note that the exact upper and lower limits on 
the luminosity are unimportant since P-yP-y that enters in Eq. ( ^ ) is peaked well between the lower and upper limits 
and is vanishingly small at the boundary. For the redshift integration we will as a default integrate between Zaihi = 
and 2max = 5, but this integration range will be, as discussed below, restricted to include the effect of resolved blazars. 
For the parameters k, rj, ^ and Af, we will use the values obtained in ||60|| , 



(46) 



where k was determined from observations of blazars that are observed both in radio and in gamma rays, rj was 
determined by requiring the number counts of blazars to be consistent with the EGRET observations and ^ and A f 
was determined from EGRET blazar observations. 



K 


= 4-10 


V 


= 1.0 


C 


= 0.03 




= 5.0 



C. Intrinsic gamma ray spectrum 



We will assume that the intrinsic gamma ray spectrum follows a power law with spectral index a, i.e. that 



E 



(47) 
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where Ef = OA is the fiducial energy at which we calculate the luminosity P^. Note that it is probably unrealistic 
to assume that the spectrum continues to be a power law to higher energies (above a few hundred GeV), instead we 
should expect a cut or a tilt in the spectrum. However, we will here for simplicity assume that there is no cut-off 
which means that we will probably overestimate the diffuse gamma ray background at high energies. 



D. Flux from a single source 

When taking resolved blazars into account we need the gamma ray flux a given blazar would produce. A blazar 
with luminosity and spectral index a at redshift z will give rise to the integrated gamma ray flux above energy 
Eth, 

Py Eth (Eth{l + z)y 1 



2'KhEfa-l\ Ef J {RoSk{r)y 

This equation is valid for Eth ^ 10 GeV since we here have neglected absorption (which is a reasonable approximation 
for low energies). With appropriate unit conversions this is the flux in units of cm^^ s^^ that should be compared 
with the EGRET or GLAST point source sensitivity. For EGRET we will use the point source sensitivity 1 • 10~^ 
cm-2 s-i i^l and for GLAST we will use 1.6 • 10"^ cm'^ s'^ [El. 



E. Distribution of spectral indices 



We have to choose a distribution function for the spectral indices, p{a). One option is to use the distribution of 
spectral indices of blazars as observed by EGRET, 



1 j:^. 



1 



(49) 



where the sum is over the observed spectral indices a; with their corresponding errors Ui. However, this is not the 
best choice of distribution function since sources with low a are easier to detect due to their harder spectrum and we 
would hence introduce a selection bias. Instead we select a distribution function of the form 



1 



(JintV 27r 



(50) 



where we fix aint and (Tint such that the predicted distribution of a for observable blazars match the observed 
distribution. The predicted a distribution as it should be observed by EGRET is given by 



Pohs{a)da 



1 47rc /•i°Sio-Pr'= 
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dlogio^^ 
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dz 
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Pintia- Aaq)p^^g{P^,z)+pintia~ Aaf)py J {P^,z) (RoSkir)) da. 



(51) 



where we only integrate over observable blazars, which is most easily done by noting that a blazar at redshift z, 
with luminosity P^ and spectral index a is observable if it would produce a flux above the EGRET point source 
sensitivity of 1 • 10"'' cm~^ s~^ integrated above 0.1 GeV. Using Eq. (^) above we can for a given P^ and a calculate 
the maximum redshift z' at which such a blazar would be observable. This would then be our upper limit for the 
z-integration, i.e. Zmin — and Zmax — z' . A^prcd in Eq. (O) is the total number of observable blazars and is given by 



prcd 



da 



dz 



1 



Hz) 



p{a - Aaq)pj^q{Pj, z) +p{a - Aaf)p^j{P-f,z) (RoSkir)) 



(52) 



We now have to choose a sample of observed blazars and fit aint and dint such that we can reproduce the observed 
distribution of a. We have followed this procedure for two samples of blazars, the first one are 27 blazars by Lin et al. 
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Spectral index, a Spectral index, a 



FIG. 10: Different a distributions, a) The solid line is the intrinsic distribution A discussed in the text, the dashed is the 
predicted EGRET observable distribution and the dash-dotted line is the observed EGRET distribution for the Lin et al. 
sample b) The same as in a) but for intrinsic a distribution B and compared with the sample in the 3rd EGRET catalog 



p3[ and the second one are the 65 blazars with determined spectral indices in the 3rd EGRET catalog |64j. Before 
we do this fit, we fix the spectral shifts of blazars in quiescent and flaring states as 



-0.05 



Aaf = 0.20 



(53) 



which are the shifts determined by Stecker & Salamon for EGRET blazars which have been observed in both 
quiescent and flaring states. For the two samples we then get 



= 2.25 
^ 0.30 



Lin et al. 



'int 



= 2.35 
= 0.30 



3rd EGRET catalog K4 



(54) 



These values are in very good agreement with the results in 1 65 . We will refer to the first and second set of parameters 
as distribution A and B respectively. In Fig. |lO|a w e compare distribution A with the predicted EGRET distribution 
and the observed EGRET distribution. In Fig. |10|b we do the same thing for distribution B. Note that both predicted 
distributions fit the two observed samples rather well, but that the sample in the 3rd EGRET catalog is shifted by 
about 0.1 compared to the Lin et al. sample. This shift is of the same order as the expected systematic uncertainty 
in the EGRET catalog. In the following, we will use distribution A as our defa ult s ince it reproduces the EGRET 
observed diffuse extragalactic background better than distribution B (see section V F below) 

The predicted number of observed blazars is given by Eq. ^ 
^pred ^42, in reasonable agreement with the observed number of 66 blazars [|64|. Note that we do not expect perfect 
agreement since we only use a simple point source sensitivity, but in reality the sensitivity is much more complicated. 
We could easily envision that it should depend on e.g. the spectral index a. Hence we are content that the agreement 
is as good as it is. Note that we also have the freedom to change the beaming parameter rj, but we choose to keep it 
fixed to 7y = 1 as given in Ref. [|60[. 



and for the two distributions we get N^^^^ = 51 and 



F. Taking resolved blazars into account 

In Eq. (^) we should only integrate over unresolved blazars. This is done in the same way as in the previous 
section, i.e. for a given luminosity Pj and spectral index a there is a given redshift z' below which the blazars will be 
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FIG. 11: The predicted diffuse gamma ray flux (multiplied by to show leatures more clearly) for EGRET, a) The predicted 
fluxes for different a distributions. Distribution C is hke A and B but with aint = 2.15. As can be seen the exact shape of 
the spectrum is fairly sensitive to the a distribution. Also shown are the EGRET measurements of the diffuse extragalactic 
gamma ray background, b) The predicted flux for different absorption models. 



resolved and above which they will be unresolved. If we let the lower limit of the redshift integration, be equal 
to this redshift z' we will only include unresolved blazars. 

In Fig. |ll| we show the predicted diffuse gamma ray fluxes for EGRET. As can be seen in a), our model reproduces 
the measured diffuse extragalactic background fairly well. To further illustrate the dependence on the a distribu- 
tion, we also show results for a hypothetical a distribution (C) with a-mt = 2.15. For this distribution, the agreement 
with the EGRET measurements is excellent (the slight difference in normalization could be fixed by slightly increasing 
the beaming parameter ry). We should note that our predictions are fairly sensitive to the exact low-a behavior of 
p(a). The higher up in energy we go, the more we sample the low-a region. In b) we show the effects of the different 
absorption models. It is clear that as soon as we go above 10-100 GeV, absorption effects are very important. Also 
keep in mind that we have not included a cut-off in the intrinsic gamma ray spectrum which would further reduce 
the fluxes at high energy. 

In Fig. [l^ we show the effect of different point source sensitivities. We see that compared to EGRET, the superior 
point source sensitivity of GLAST will reduce the diffuse gamma ray background with roughly a factor of two. Note 
however, that the angular resolution of GLAST will make the point source sensitivity worse at lower energies (or 
rather, larger spectral indices a), an effect we have not included here. We expect that this effect would make the 
predicted background for GLAST slightly higher at low energies than shown in the figure. In Fig. p^ b we show the 
effect of different absorption models for the predicted GLAST gamma ray background. 



G. Uncertainties 



In this section we have produced a derivation of the expected diffuse gamma ray background assuming that it is 
due to unresolved blazars. There are many uncertainties involved. First of all, it is not known whether blazars are the 
only sources relevant to compute the background. The energy spectrum of the blazars is also not very well known, 
i.e. there could be a cut-off at high energies (and even if the spectrum is a power law, the distribution of spectral 
indices is uncertain). Even the luminosity function is rather uncertain and the assumption of the relation between 
the gamma and radio luminosity functions cannot be tested until the sample of blazars measured in both gamma 
and radio increases. The parameters of the model we discussed are also quite uncertain, and, as already mentioned, 
gamma ray absorption introduces further uncertainties, especially at high energies. In spite of all these uncertainties, 
the agreement we find between our prediction and EGRET data is quite good, and gives some credibility to our 
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FIG. 12: The predicted diffuse gamma ray flux (multiplied by to show features more clearly) for GLAST. a) The predicted 
fluxes compared to the EGRET measurements, b) The predicted fluxes for different absorption models. 



estimate of the background for GLAST at higher energies. We have chosen as our default model the a distribution A, 
which reproduces both the measured a distribution and the EGRET energy spectrum satisfactory, and the absorption 
model of Primack et al. Keep in mind though that, above ~ 100 GeV, the uncertainties can be as large as a factor of 
a few. 



VI. APPLICATIONS TO SUPERSYMMETRIC DARK MATTER 



A. A few examples in a specific particle physics setup 

So far, we have kept the discussion as general as possible, without specifying the exact identity of the WIMP 
making up the dark matter. To gauge the possibility of detecting a gamma-ray signal in a realistic scenario, we now 
investigate one of the prime candidates for dark matter: the lightest supersymmetric particle (LSP) in the MSSM - 
the minimal supersymmetric extension of the Standard Model of particle physics. If R-parity is conserved, the LSP is 
stable; furthermore, its coupling with lighter standard model particles ensures that a population of such particles is 
present in the early Universe, with its density set by thermal equilibrium. The freeze out from equilibrium is roughly 
set by the LSP thermally averaged annihilation cross section; as sketched in Eq. (p9|), a weak interaction strength 
coupling ensures that WIMPs have a thermal relic abundance of the order of the critical density: this is naturally the 
case if the LSP has zero electric and color charges. 

We thus take as our template WIMP dark matter candidate the lightest neutralino, Xi, in the MSSM (see |]l| for 
a recent review). Xi is defined as the lightest mass eigenstate obtained from the superposition of four spin-1/2 fields, 
the Bino and Wino gauge fields, B and W^, and two neutral CP-even Higgsinos, and i/j • 

X? = NnB + Ni2W^ + NisHi + N^H^ . (55) 

There are large regions in the MSSM parameter space where Bino-like LSPs or neutralinos with relevant Higgsino 
components have a thermal relic abundance of the right order to account for the dark matter, see, e.g., [|66| . We 
have used the DarkSUSY package |^ to scan extensively the parameter space and generate a large archive of such 
models. We select those models that do not violate present accelerator and astrophysical limits and study what 
is the typical gamma-ray yield, both for the continuum and monochromatic spectra (in the MSSM there are two 
line signals allowed: 77 and Z7). With DarkSUSY |^ we calculate the relic density by numerically solving the 
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FIG. 13: Extragalactic gamma-ray flux (multiplied by E'^) for two sample thermal relic neutralinos in the MSSM (dotted 
curves), summed to the blazar background expected for GLAST (dashed curve). Normalizations for the signals are computed 
assuming halos are modelled by the Moore profile, with 5% of their mass in substructures with concentration parameters 4 
times larger than c^ir as estimated with the Bullock et al. toy model. 



Boltzmann equation properly taking resonances, thresholds and coannihilations (between the lightest neutralinos and 
other neutralinos and charginos) into account []66|. 

We focus first on neutralinos with a relic abundance fl^h^ in the interval 0.1 to 0.2, corresponding to our preferred 
cosmology ^Im ~ 0.3 and h ^ 0.7. There are models with (TV27 ^ 10~^^cm'^s~"'^ over the whole mass range from 50 
GeV up to a few TeV. We consider two sample cases and plot the corresponding extragalactic 7-ray spectra in Fig. |l3 
(dotted lines). The first model has M-^ = 76 GeV, relatively low total annihilation cross section crv = 6.1 • 10~^^cm^s 
but large branching ratios into photon states, 627 = 6.1 • 10~^ and bzj = 5.2 • 10~^. The other has = 171 GeV, 
larger annihilation rate av — 4.5 • 10~^^cm'^s~^ but 627 = 5.2 • 10~* and negligible branching ratio into the Zj final 
state. The normalization of the flux is set by assuming that dark matter structures are described by the Moore profile, 
with concentration parameters as computed with the Bullock toy model, and assuming the presence of a moderate 
amount of substructure, / = 5%, with a factor of 4 enhancement in Cyir- Under these assumptions, the neutralino 
induced 7-ray flux is at the level of the diffuse background from unresolved blazars (a distribution A) expected in 
GLAST, with the peak from the monochromatic emission significantly above it (dashed curves refer to the background 
only, solid curves to the sum of signal plus background). 

The condition il-^h^ > 0.1 sets an upper bound on to the total annihilation cross section and hence, indirectly, 
an upper bound on the strength of the monochromatic channels; such states however are not the dominant modes 
and therefore a lower bound does not follow from imposing fl^h^ < 0.2: there are cases where the Xi is compatible 
with being a good dark matter candidate, but the monochromatic flux is negligible. An opposite conclusion holds for 
the continuum components: there are cases in which the gamma-ray yield can be slightly larger than the one for the 
= 171 GeV model, but very small regions in parameter space where the yield is significantly smaller than for the 
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FIG. 14: Extragalactic gamma-ray flux (multiplied by E'^) for two sample non-thermal dark matter candidates arising in the 
AMBS scenario (dotted lines) compared to the expected background (dashed curve). Annihilations cross sections are in these 
cases larger than for the models displayed in Fig. Iisl however a different normahzation for the fluxes is implemented here: we 
consider the case for halos modelled by the NFW profile, no substructures and concentration parameters inferred from the 
BuUock et al. toy model. 



model with = 76 GeV we show. 

If we remove the constraint on il^ the picture can change drastically. In particular, there are several schemes in 
which the LSP relic abundance today is not set by its thermal relic density. One example is the case for Wino or 
Higgsino-like neutralinos in the version of the MSSM with anomaly mediation for supersymmetry breaking (AMSB). 
This scheme induces a mapping into regions in the MSSM parameter space in which the thermal relic abundance is 
negligible; however, an additional "non-thermal" relic source is present due to decays into neutralinos of gravitinos or 
moduli fields, fields that parameterize a flat direction of the theory and that dominate the energy density in the early 
Universe One can show that, in this context, the total annihilation rate, as well as cross sections in the 2"f 

and Zj final states, are forced to be very large ||7^. Two examples are shown in Fig. one model has = 92 GeV, 
av = 2.5 • 10"24(,j^3g-i^ ^ 12 ■ 10-3 and bzy = 2.2 • lO^^; the second = 180 GeV, crv = 2.2 • 10"24j,j^3g-i^ 
627 = 1.8 • 10~'^ and bz^y = 5.1 • 10-^. The normalization of the two extragalactic 7-ray fluxes is set assuming NFW 
halo profiles with no substructure and concentration parameters as computed with the Bullock toy model. Had we 
chosen the Moore profile rather than NFW, the predicted fiuxes would be one order of magnitude larger, hardly 
compatible with the extragalactic flux as inferred from EGRET data. Note that a fiux at roughly the same level is 
expected implementing the Burkert profile, hence the detectability of this signal is not linked to having a singular 
halo profile describing dark matter halos. 
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FIG. 15: Approximate 3cr sensitivity curves for the GLAST telescope to search for a component in the extragalactic gamma- 
ray flux induced by WIMP annihilations into monochromatic photons. The sensitivity curves are plotted in the plane WIMP 
mass (coinciding with the energy peak in the induced flux) versus twice the annihilation rate into two photons, and for four 
configurations to estimate the normalization of the flux (the highest and lowest dashed curves correspond, respectively, to the 
choice in Fig. |l^ and Fig. ^3|). Also shown is the range of the predictions of VfJ2~i for neutralinos in the AMSB scenario, and the 
upper limit to it in case of thermal relic neutralinos in the MSSM, assuming their relic abundance is either in the cosmologically 
preferred mass range 0.1 < 0.^1? < 0.2, or in the less restrictive range often considered 0.025 < f2^/i^ < 0.2. 

B. Sensitivity in upcoming measurements 

It is not straightforward to estimate the smallest WIMP induced component GLAST will be able to disentangle 
from the background. A firm statement about the possibility to single out the yield with continuum energy spectrum 
will be possible only when higher precision measurements will allow to characterize better the level and the spectral 
features of the background. The monochromatic component has a much better signature and might be unambiguously 
identified. We make an attempt to make a rough estimate of the sensitivity curves for a GLAST-like instrument under 
a few schematic assumptions. 

We assume the instrument has a peak effective area of 8000 cm'^ at energies above 10 GeV and an average energy 
resolution of 15% |5^. We take an exposure of 4 years, mapping the whole sky except for the regions already excluded 
in the EGRET analysis Q, i.e. the galactic plane |6| < 10°, and the bulge |/| < 40° and |6| < 30°, with an average 
effective area which is about 20% of the peak area. We set up a procedure to check if we can discriminate the 
spectrum of a line signal superimposed on the background from the spectrum of the background only. The analysis is 
performed assuming a given normalization for the WIMP flux and keeping as free parameters the value of the WIMP 
mass and annihilation cross section uvi-^ . For each parameter choice, we sum to this flux the diffuse background from 
unresolved blazars (a distribution A) with the normalization computed in the previous Section and already shown in 
Figs. O-ll^. We then perform a binning of the spectrum above 10 GeV, optimizing the bin width as a function of 
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the number of events in each bin and checking that we have at least 10 events per bin. Naively, the statistical error 
in each bin would be the square root of the number of events in the bin; at second thought though, the extragalactic 
background component will be obtained after subtracting point sources and the diffuse galactic emission, with a non 
trivial propagation of errors we cannot easily retrace here. We make a simplifying assumption at this stage, expecting 
just a rough estimate of the true sensitivity curves. Suppose the main component one has to fight against is due to 
diffuse galactic 7-rays; such a component can be removed after assuming a model for diffuse emission and should be, 
on average, about an order of magnitude larger than the extragalactic component [ p3[ . We mimic this subtraction by 
assuming that the error in each energy bin is the square root of the number of events in the bin multiplied by 10. We 
then use the criterion to discriminate whether or not the obtained distribution of events can be fitted at 3 ct with 
a background component only with fixed spectral shape but arbitrary normalization. 

The corresponding sensitivity curves are shown in Fig. ^ in the plane neutralino mass - twice the 27 annihilation 
rate. Each curve corresponds to a different normalization for the extragalactic flux: from the bottom up, case for Moore 
profile halos with substructure introduced already in Fig. ^ the case for Moore profile halos with no substructure 
and concentration parameters as computed with the Bullock et al. toy-model or with the ENS model, and, finally, 
the case for NFW halos with no substructure and Cvir as in the Bullock et al. model. Also shown in the picture are 
the span in the predictions for av2-f in the AMSB scenario i fTOf (green lines marks the upper and lower limits for a 
given mass), and approximate upper limits in the case of MSSM thermal relic neutralinos with relic abundance in the 
preferred range 0.1 < U^h^ < 0.2, or in the less restrictive range often considered 0.025 < fi^/i^ < 0.2, as deduced 
from our sampling of the parameter space. As can be seen, depending on the configuration one considers, there is 
a fair chance that the monochromatic 7-ray flux will be disentangled from GLAST data. The four models we have 
considered in Figs. |l3| and |lj lie all above the corresponding sensitivity curves. 

The same sensitivity curve can be applied to the case of the line signals generated in the Z7 channel by replacing 
M-^ on the horizontal axis with the energy of the peak E — M^{1 — Aff /4M^) and assuming the quantity on the 

vertical axis is avz-y ■ 



C. Comparison with other signals 

It is not straightforward to compare the dark matter signal we have presented here with other indirect signals 
which were proposed soon after the idea of WIMP dark matter was raised, two decades ago. Most analyses have been 
devoted to the study of the detectability of gamma-rays produced in the halo of our own Galaxy or of antimatter 
generated by WIMP pair annihilations taking place in our local environment (say within a few kpc, the exact number 
depends on the model for propagation of charged cosmic rays), refining the original proposals, see ||l| for a detailed 
reference list. 

As already mentioned, the prospects for detecting gamma-rays produced in the Milky Way are much more tightened 
to assumptions on the distribution of dark matter WIMPs in the galactic halo. The monochromatic flux generated 
by the sample MSSM models displayed in Fig. |l^ and Fig. |l^ in the Galactic center region is within the sensitivity 
of GLAST or the upcoming generation of ground based air Cherenkov telescopes (see, e.g., the analysis in Ref. [ [70| ) 
if indeed the dark matter density profile is singular, respectively, as in the Moore et al. or the NFW halo all the way 
to the central black hole (or maybe even steeper than that, see the possible enhancement induced by the black hole 
formation discussed inj[7l|, but take into account also the opposite conclusions drawn in, e.g., Refs. ^). As 
suggested also by Fig. ^, even a slight depletion in the central density can change drastically this conclusion. 

We checked also that the four sample MSSM models we introduced, with the halo profiles considered in the 
corresponding figures, do not generate a continuum spectrum component which exceeds the flux measured by the 
EGRET telescope [Q . A comparison with the Galactic flux at high latitudes in a configuration with clumps in the 
halo is much more uncertain. The flux is dominated by eventual nearby clumps, depending critically on the actual 
physical realization which happens to correspond to the Milky Way: we recall, on the other hand, that the signal 
we propose is obtained as the sum of many unresolved sources, i.e. we are automatically making an average over a 
set of possible configurations. The chance for GLAST to resolve single clumps, for the sample configurations with 
clumps considered in Fig. |l^, can be read out of Fig. ||. The models with neutralino masses = 76, 171 GeV, have, 
respectively, A'^"'^ = 27.4 and 39.9; hence dotted lines corresponding to the Moore profile in Fig. ^ should be rescaled 
along the vertical axis by, respectively, a factor of 0.096 and 0.45. 

Limits from charged cosmic ray data are also model dependent, as again the dark matter signal is dominated by 
local sources; dark matter candidates may be excluded in some configurations, but allowed in others. Notice also 
that, especially if one focusses on the monochromatic gamma-ray component, such a signal is very weakly correlated 
to the production rate of e.g. antiprotons and positrons, see, e.g., ||75|. 
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VII. CONCLUSIONS 



We have studied predictions and the observabihty of the diffuse gamma ray signal from WIMP pair annihilations in 
external halos. We have found configurations that imply signals at a detectable level for GLAST, the upcoming gamma- 
ray space telescope, both for non-thermal dark matter neutralinos, such as in the anomaly-mediated supersymmetry 
breaking model, and for thermal relic neutralinos in the MSSM. The key ingredient to show that detectable fluxes may 
arise, which was neglected in early estimates, is the picture, inspired by the current theory for structure formation and 
by N-body simulation results, that dark matter clusters hierarchically in larger and larger halos, with light structures 
more concentrated than more massive ones. For dark matter candidates in the AMSB scenario our conclusion holds 
independently of further assumptions on the dark matter distribution inside halos. Pair annihilation cross sections 
for thermal relic WIMPs are generally smaller; this however can be compensated by the enhancement in the flux one 
finds if, as suggested by results of simulations, we assume that dark matter profiles are singular and contain small 
dense substructures. 

If the branching ratio for WIMP annihilation into monochromatic gamma rays is significant (about a few times 10~^ 
or larger), the induced extragalactic flux shows a very distinctive feature, the asymmetric distortion of the line due to 
the cosmological redshift and its sudden drop at the value of the WIMP mass. The component with continuum energy 
spectrum can be at the level of background components but has less distinctive features: the flux is characterized by 
the "7r° bump", rather than by a spectral index, with the peak shifted to energies lower than M^/2 and the width set 
by the WIMP mass. Once a better measurement of the background will be available, it will be possible to address 
the question of whether or not this signal can be disentangled from other eventual components. 

We have discussed in detail how our predictions depend on assumptions on the Cosmological model and the structure 
formation picture applied. Unless one introduces drastic changes, such as a large warm dark matter component, the 
cosmological parameters in the CDM setup do not play a major role; results are mainly sensitive to erg with about a 
factor of two uncertainty. Larger indeterminations, of the order of a factor of a few, are introduced when estimating 
the scaling of the concentration parameter with halo mass, as extrapolations with toy models out of the mass range 
of N-body simulation results are needed. The functional form of the dark matter profile in single halos introduces a 
factor of 10 uncertainty, going from the case of a cusp in the Moore profile to the case of non-singular profiles; 

that uncertainty is much smaller than, e.g., the one induced on the estimate of the flux from the center of our own 
Galaxy. The presence of substructures inside halos may provide a factor of a few enhancement in the flux, but this 
effect is more difficult to address: we have presented a simple and rather generic setup, which will be possible to refine 
when further information on halos will be provided by higher resolution numerical simulations. 

Issues related to the estimate of the background are very important as well. We have presented here a novel estimate 
of the expected background from unresolved blazars in GLAST, exploiting recent data and discussing critically the 
uncertainties involved, including the role played by gamma absorption. 

Concluding, the present analysis has been devoted to examining in detail an idea that three of the authors have 
recently presented in a short letter This work provides further support for such proposal, with exciting perspectives 
for upcoming measurements. 
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Swedish Research Council for support. 
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